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TECHNICAL  FINAL  REPORT 
April  1,  1979  -  September  30,  1981 


The  present  report  summarizes  the  research  effort  by  the 
Materials  Sciences  Corporation  on  the  analysis  of  temperature- 
dependent  stress-strain  relations  of  metal-matrix  composites 
(MMC)  during  the  period  April  1,  1979  through  September  30, 

1981.  The  work  done  is  part  of  a  comprehensive  program  to  pro¬ 
vide  analytical  tools  for  the  evaluation  of  thermomechanical  pro¬ 
perties  and  internal  stresses  of  unidirectional  MMC.  This  pro¬ 
gram  consisted  of  the  following  parts: 

1.  Analytical  determination  of  one-dimensional  temperature 
dependent  stress-strain  relations  on  the  basis  of  micro¬ 
mechanics  . 

2.  Determination  of  temperature-dependent  stress-strain 
relations  for  combined  states  of  stress  in  terms  of 
one-dimensional  stress-strain  relations  on  the  basis 
of  macro-incremental  theory. 
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INTRODUCTION 


The  usual  unidirectional  metal  matrix  composite  (MMC)  con¬ 
sists  of  an  aluminum  alloy  matrix  and  aligned  carbon,  graphite 
or  boron  fibers.  The  metal  matrix  is  a  temperature  dependent 
elasto-plastic  isotropic  material.  The  fibers  are  anisotropic 
elastic  throughout  the  temperature  range  but  their  properties 
may  be  temperature  dependent.  They  may  be  described  as  trans¬ 
versely  isotropic  elastic  with  respect  to  their  longitudinal 
axes. 

Metal  matrix  composites  are  considerably  more  expensive 
than  polymer  matrix  composites.  The  main  reasons  for  their  de¬ 
velopment  are  their  better  retention  of  strength  and  stiffness 
at  elevated  temperatures.  For  purposes  of  efficient  designs 
with  such  composites,  their  mechanical  behavior  must  be  quan¬ 
titatively  known.  The  composite  is  evidently  an  anisotropic 
temperature  dependent  elasto-plastic  material.  Consequently, 
all  applied  stresses  and  temperatures  will  produce  interactive 
non-superposable  deformations.  In  addition,  the  composite  is 
history  dependent.  This  implies  that  different  loading  histories 
resulting  in  the  same  end  state  (e.g.  axial  load-torsion-heatina , 
torsion-axial  load-heating,  etc.)  will  produce  different  deforma¬ 
tions.  It  follows  that  an  experimental  program  to  define  the 
stress-strain  relations  would  require  an  enormous  number  of  ex¬ 
periments  and  the  interpretation  of  the  results  would  be  ex¬ 
tremely  difficult. 

Another  alternative  is  to  devise  an  analytical  treatment 
for  determination  of  the  stress-strain  relations.  This,  however, 
also  involves  severe  difficulties.  Analysis  of  any  model  of  a 
fiber  composite  can  only  be  carried  out  numerically  and  the 
amounts  of  computer  time  required  make  such  procedures  quite 
expensive.  It  is  in  principle  possible,  as  will  be  shown  in 
this  work,  to  carry  out  such  numerical  analyses  for  any  load- 
temperature  program  but  here  there  arises  the  same  difficulty 


as  with  the  experimental  program.  The  number  of  variables  is 
large  and  load-temperature  histories  are  of  infinite  variety. 

Thus,  purely  numerical  analysis  is  as  expensive  and  as  futile  as 
purely  experimental  investigation. 

A  possible  solution  to  these  difficulties  is  to  restrict 
numerical  analysis  to  simple  loading  cases  in  which  only  one  stress 
component  is  applied,  e.g.  simple  tension  or  pure  shear,  and 
to  establish  a  global  theory  for  the  stress-strain  relations  for 
complicated  loadings  which  incorporates  the  information  contained 
in  the  one  dimensional  stress-strain  relations,  obtained  numeri¬ 
cally.  (Such  a  global  theory  can  also  be  based  on  experimentally 
obtained  one  dimensional  stress-strain  relations.)  This  method 
will  be  followed  in  the  present  work. 

The  problem  of  the  analytical  determination  of  elasto-plastic 
stress-strain  relations  has  received  repeated  attention  in  the 
literature.  A  general  qualitative  discussion  of  elasto-plastic 
composites  has  been  given  in  reference  1.  Limit  analysis  treat¬ 
ments  which  are  only  concerned  with  prediction  of  ultimate  loads 
for  ideally  plastic  matrix  have  been  given  in  references  2  and  3. 

A  numerical  treatment  to  compute  the  stress-strain  relation  for 
transverse  uniaxial  stress  for  a  limited  strain  range  has  been 
given  in  reference  4  based  on  the  idealized  geometry  of  a  square 
array  fiber  model.  A  similar  approach  but  including  interaction 
between  shear  and  transverse  load  has  been  given  in  reference  3. 

A  micromechanics  stress  analysis  for  a  square  array  has  been 
given  in  reference  6.  The  self-consistent  scheme  approximation 
has  been  utilized  in  reference  7  for  the  case  of  rigid  fibers 
and  elasto-plastic  matrix  characterized  by  isotropic  deforma¬ 
tion  theory.  This  was  generalized  in  reference  8  to  include 
fiber  elasticity;  all  this  for  transverse  loading  only.  Signifi- 
can  work  to  predict  initial  yield  surfaces  has  been  given  in 
references  9,  10,  11  and  12,  including  temperature  change  but 
primarily  for  axisymmetric  states.  The  general  form  of  the 
elasto-plastic  isothermal  stress-strain  relations  of  a  unidirec- 


tional  fiber  composite  consisting  of  elastic  isotropic  fibers  and 
ideally  plastic  matrix  has  been  discussed  in  reference  13. 

All  of  the  treatments  available  describe  important  special 
cases  but  they  are  not  sufficient  to  define  the  needed  tempera¬ 
ture  dependent  stress-strain  relations  of  a  unidirectional  MMC, 
neither  for  the  one  dimensional  nor  for  the  three  dimensional 
case . 

The  following  investigation  is  divided  into  three  main  parts. 

In  the  part  METHOD  OF  ANALYSIS:  ONE  DIMENSIONAL  EFFECTIVE  STRESS- 
STRAIN  RELATIONS,  we  formulate  the  elasto-plastic  problem  of  the 
hexagonal  array  of  fibers  in  terms  of  two  kinds  of  two  dimensional 
boundary  value  problems  for  typical  repeating  elements.  This  for¬ 
mulation  is  then  used  to  numerically  obtain  effective  temperature 
dependent  stress-strain  relations  for  one  dimensional  loadings 
such  as  transverse  tension,  transverse  shear  and  axial  shear  by 
means  of  the  ANSYS  computer  code. 

In  the  part  THREE  DIMENSIONAL  STRESS-STRAIN  RELATIONS. 
ISOTHERMAL  CASE,  we  use  global  plasticity  considerations  to  con¬ 
struct  isothermal  incremental  stress-strain  relations  of  a  trans¬ 
versely  isotropic  unidirectional  fiber  composite  in  terms  of  known 
one  dimensional  stress-strain  relations. 

Finally,  in  the  part  TEMPERATURE  DEPENDENT  STRESS-STRAIN  RE¬ 
LATIONS,  we  extend  these  results  to  temperature  changes  resulting 
in  a  set  of  effective 'stress-strain  relations  incremental  in  stress, 
strain  and  temperature. 


-3- 


METHOD  OF  ANALYSIS:  ONE  DIMENSIONAL 
EFFECTIVE  STRESS-STRAIN  RELATIONS 


When  a  homogeneous  body  consisting  of  elasto-plastic 
material  is  subjected  to  the  traction  rate  boundary  conditions 

T°(S)  =  o  ?  j  n  .  (2.1) 

where  i'L  are  space  constant  stress  rates  and  n_.  are  components 
of  the  outward  unit  normal,  the  stress  rates  throughout  are  the 
constant  "•  °  j .  If  (2',  1)  are  applied  to  the  surface  of  an  elasto- 
plastic  heterogeneous  body,  such  as  the  fiber  composite  under 
consideration  here,  it  follows  from  the  average  stress  rate 
theorem  (ref.  14)  that  average  stress  rates  are  given  by 


If  the  internal  phase  geometry  of  the  component  is  statistical ly 
homogeneous,  the  stress  and  strain  fields  associated  with  (2.1) 
are  also  statistically  homogeneous.  Then  (2.2)  can  be  inter¬ 
preted  in  terms  of  averages  over  a  representative  volume  element 
(RVE)  . 

Now  let  the  body  be  subjected  to  the  constant  boundary 
temperature  rate 


b  (S)  =  bo  (S)  . 


(2.3) 


It  follows  trivially  from  the  equations  of  steady  state  un¬ 
coupled  heat  conduction  that  the  temperature  rate  field  through¬ 
out  the  heterogeneous  body  is  the  constant  : 0  and  thus  the 
temperature  is  space  constant  at  all  times. 


-4- 


In  a  temperature  dependent  elasto-plastic  material  the 
relation  between  strain  rates,  stress  rates  and  temperature 
rate  is  linear.  It  can  be  shown  that  if  the  constituents  of 
a  composite  have  this  property,  the  same  is  true  for  the  re¬ 
lation  of  average  strain  rates  to  average  stress  rates  and 
temperature  rate  when  the  composite  is  subjected  to  (2.1) 
and  (2.3).  Therefore, 


C*  +  * 

ijkl  kl  i  j  '  ° 


(2.4) 


where  C?.,,  and  are  complicated  functions  of  average  stress, 

1 3  k  i  i  3 

strain  and  temperature  history  and  of  the  internal  phase  geometry, 
but  not  of  average  stress  rates,  strain  rates  or  temperature 
rate.  The  relations  (2.4)  are  the  effective  stress-strain  re¬ 


lations  of  the  composite.  Their  fo^m  will  be  discussed  later 
in  this  work.  At  present,  attention  will  be  directed  towards 
stress-strain  response  for  one  dimensional  temperature  depend¬ 
ent  loadings. 

We  now  consider  the  specific  case  of  a  unidirectional  fiber 
composite  in  the  form  of  a  long  cylindrical  specimen  (figure  1). 
The  matrix  is  temperature  dependent,  isotropic  and  elasto-plastic, 
and  the  fibers  are  elastic  transversely  isotropic.  This  is 
an  appropriate  description  of  aluminum  matrix  containing  carbon 
or  graphite  fibers.  To  compute  effective  stress-strain  relations, 
it  is  necessary  to  find  the  average  strain  rates  produced  by  the 
average  stress  rates.  This  requires  analysis  of  strain  and  stress 
rate  fields  ‘subject  to  displacement  and  traction  rate  continuity 
at  fiber/matrix  interfaces,  the  analysis  being  carried  out  in¬ 
crementally  for  specified  external  loading  histories.  Since 
the  matr-ix  is  temperature  dependent  elasto-plastic,  such  analy¬ 
sis  for  any  given  model  of  a  fiber  composite  can  only  be  per¬ 
formed  numerically. 
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All  of  the  problems  to  be  considered  here  can,  fortunately, 
be  formulated  in  two  dimensions.  For  this  purpose,  let  the 
strain  rates  (2.2)  be  split  as  follows 


'a  =  '.P  +  '.a 

"ij  ij  “ij 


(2.3) 


(a) 


(b) 


(2.6) 


It  may  be  shown  that  when  r”.  in  (2.1)  have  the  form  (2.6a),  the 
problem  of  the  determination  of  the  stress  and  strain  rates  in 
fibers  and  matrix  can  be  formulated  in  terms  of  generalized  plane 
strain  which  implies  that 


L11 

const 

(a) 

1  *  (.*' 

c„,  <x2'X3) 

(b) 

>,-■=2,3 

(2.7 

'12 

C  1 3  =  0 

(c) 
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in  fibers  and  matrix.  Such  problems  can  be  solved  in  terms 
of  two  dimensional  stress  rate  fields  with 


12 


(2.8) 


When  c  V  j  have  the  form  (2.6b),  the  problem  can  De 
formulated  in  terms  of  anti-plane  strain  which  implies  that 


c12  i12(x2,x3> 

C13  =  c13(x2'x3’ 


(2.9) 


and  all  other  strain  rates  vanish.  Such  problems  can  be  solved 
in  terms  of  a  two  dimensional  stress  rate  field  for  which 


’ll  "  '22  :33  ~  c23  0 


(2.10) 


and  thus  only  the  two  shear  rates  c ^  and  '^3  remain,  which 
are  functions  of  *2,x3  only. 

One  dimensional  loadings  of  interest  are  uniaxial  average 
stress  in  the  fiber  direction,  transverse  uniaxial  stress,  trans¬ 
verse  shear  and  axial  shear.  The  first  three  all  fall  into  cate¬ 
gory  (2.6a)  while  the  last  is  in  category  (2.6b).  For  stress  in 
the  fiber  direction,  when  the  fibers  are  high  modulus  carbon, 
graphite  or  boron,  it  may  be  safely  assumed  that  the  axial  av¬ 
erage  strain  is  elastic  since  in  this  case  fiber  and  matrix 

Ejj  strains  must  be  equal  and  the  stiff  fibers  inhibit  plastic 
matrix  flow.  In  the  other  cases,  however,  substantial  average 
plastic  strains  can  develop  as  a  consequence  of  matrix  plasticity. 


For  the  purpose  of  computation  of  effective  stress-strain 
relations,  it  is  necessary  to  describe  the  fiber  composite 
by  a  suitable  geometrical  model.  There  exist,  at  the  present 
time,  only  two  kinds  of  models  which  permit  exact  analysis:  (1) 
the  composite  cylinder  assemblage  (CCA) (refs.  14,15)  shown  in 
figures  2,  3  which  has  been  employed  successfully  to  find 
effective  properties  for  various  kinds  of  linear  physical 
behavior,  in  closed  analytical  form,  and  (2)  periodic  arrays  of 
identical  circular  fibers  which  have  been  analyzed  numerically. 
Unfortunately,  the  CCA  is  only  of  limited  usefulness  for  the 
present  case  of  elasto-plastic  temperature  dependent  matrix 
since  it  does  not  lend  itself  to  evaluation  of  transverse  nor¬ 
mal  or  shear  stressing  while  analysis  of  axial  shearing  would 
require  complicated  numerical  procedures.  Since  the  fiber 
composite  is  transversely  isotropic  in  the  macrosense,  we 
shall  choose  as  our  model  a  hexagonal  array  of  identical 
circular  fibers.  This  model  fulfils  the  necessary  initial 
transverse  isotropy  requirements.  Note  that  the  model  of  a 
square  array  of  circular  fibers,  which  has  often  been  used  for 
modelling  fiber  reinforced  materials,  is  not  transversely  iso¬ 
tropic  and  is,  therefore,  not  suitable.  Elastic  analyses  have 
shown  that  the  CCA  model  and  the  hexagonal  array  give  results 
which  are  numerically  very  close.  This  leads  to  the  specula¬ 
tion  that  the  results  have  wider  applicability  than  the  re¬ 
stricted  nature  of  the  model  implies  and  are  probably  valid 
for  any  transversely  isotropic  arrangement  of  separated  cir¬ 
cular  fibers. 

The  hexagonal  array  is  shown  in  figure  4.  A  typical  rec¬ 
tangular  repeating  element  is  indicated  in  the  figure.  The 
whole  array  can  be  built  up  of  such  elements  and  the  states 
of  stress  and  strain  in  these  elements  are  identical.  Actually 


the  rectangular  element  can  be  separated  into  two  repeating 
triangular  elements  as  shown  in  the  figure,  but  for  the  pur¬ 
pose  of  analysis  the  rectangular  element  is  more  convenient. 


GENERALIZED  PLANE  STRAIN 

We  first  consider  case  (2.6a)  for  some  constant  tempera¬ 
ture.  Without  loss  of  generality,  this  loading  can  be  written 


(2.11) 


where  O3  are  the  principal  transverse  stress  rates.  Be¬ 

cause  of  the  transverse  isotropy  of  the  model,  the  average 
strain  rates  produced  by  (2.11)  do  not  depend  on  the  direction 
of  transverse  axes  and  thus  the  principal  axes  can  be  taken 
as  the  Cartesian  reference  system. 

It  follows  from  symmetry  that  the  deformed  shape  of  the 
repeating  rectangle  remains  a  rectangle.  Without  loss  of 
generality,  the  normal  displacements  of  two  adjoining  sides 
may  be  taken  as  zero.  Consequently,  the  rectangle  is  subject 
to  the  following  boundary  conditions  on  its  contour  C 


°23 (C> 


u2(0,x3)  =  0 


u 2  (a2,x3)  =ri2a2 


(2.12) 


u3  (x2 , 0)  =  0  U3  (x2,a3)  =ti  a 
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A 


The  constants  n2  t-j  have  a  simple  and  important  physical 
interpretation.  Application  of  the  average  strain  rate  theorem 
(ref.  14)  in  the  two  dimensional  version  yields 


-  1  f  • 

13  h  J  13 


+  u  .n  .  )  ds 
3  1 


(2.13) 


where  A  is  the  area  of  the  rectangle.  Consequently, 


(2.14) 


n3  '  £  3  3 


In  order  to  define  the  boundary  value  for  the  rectangle, 
the  parameters  and  must  be  expressed  in  terms  of  the 
stress  rates  (2.11).  To  do  this,  we  exploit  the  rate  linearity  of 
rate  independent  plasticity.  Suppose  that  the  current  states  of 
displacement,  strain  and  stress  in  the  repeating  rectangle  are 
u,  e  and  £.  Now  consider  the  two  boundary  value  problems  shown 
in  figures  5a  and  5b  where  each  rectangle  is  in  the  same  current 
state  prior  to  displacement  increment.  If  the  solutions  to 
these  boundary  value  problems  are  identified  by  superscripts 
I  and  II,  it  follows  from  rate  linearity  that 


U  =  f^U1  + 


■  I  .  'II 

r'2l  n  3-1 


(2.15) 


-  „  *1  .  'II 

-  n2o  +  n3o 


Define  average  stress  rates 


-  =  1  fC 

0 22  a  /  u22  dx„ 

3  0  J 

°33  =  tin 


2. 

O33  dx2  • 


(2.16) 


Since  the  states  of  stress/stress  rate  are  the  same  in  all 
repeating  elements  except  for  a  boundary  layer  near  the  external 
boundary,  it  follows  from  equilibrium  that  the  stress  rates  (2.16) 
are  equal  to  (2.11).  Then  by  averaging  (2.15c) 


n2J22 


-II  _  'o 
"  3't22  °2 


(2.17) 


This  defines  the  constants  <  r, 3,  i.e.  the  average  strain 
rates  in  terms  of  stress  averages  obtained  from  the  two  bound¬ 
ary  value  problems  and  thus  the  solution  (2.15)  for  the  repeating 
rectangle  is  determined.  Note  that  r  r 3  are  functions  of  current 
state  variables.  Also  note  that  (2.15-17)  can  be  expressed  in 
terms  of  increments  instead  of  rates  and  would  be  used  in  this 
form  in  a  numerical  procedure.  Once  the  increments  du,  d_t  and 
da  are  known,  they  are  added  to  the  current  state  and  the  same 
solution  procedure  can  be  carried  out  again  for  new  increments 
daj  and  da^. 

Note  that  because  of  the  generalized  plane  strain  formulation 
the  solution  includes  the  unknown  constant  strain  rate  £.3  (2.7a). 
This  is  found  from  the  requirement  that 


•  -o  !  ra2fa3. 

°11  J1  A  J  J  °11 

0  0 


dx2dx3‘ 


(2.18) 


The  average  shear  rates  and  "3^  vanish  identically  since 

0^2  =  =  0  in  the  generalized  plane  strain  solution. 

This  procedure  defines  relations  among  the  average  stress 


rates  a  22’  ^33  and  avera9e  strain  rates  722  and  e33  at  cur¬ 

rent  states.  It  is  thus,  in  principle,  possible  to  obtain  the 
average  strains  for  specified  average  loading  history  in  incre¬ 
mental  fashion.  If  in  particular  o33  =  o°3  =  0,  then  the  pro¬ 
cedure  will  define  the  effective  stress-strain  relation  for 
uniaxial  transverse  stressing  in  the  X2  direction,  which  by  the 
transverse  isotropy  is  also  the  stress-strain  relation  for 
uniaxial  stressing  in  the  x3  direction.  The  effective  trans¬ 
verse  tangent  modulus  E*t  is  defined  by 


ETt(' 


22 


)  = 


22 


22 


22 


To  find  the  stress-strain  relation  in  transverse  shear, 
it  is  convenient  to  use  the  special  case  of  (2.11) 


°2  _  ~°3 


0°  =  0 


This  results  in  average  strain  rates 


L22  n2  "  L 33  ' 

An  element  rotated  by  45°  with  respect  to  the  x2  axes  will  be  in 
a  state  of  average  pure  shear  stress  i°  and  shear  strain  VT  de¬ 
fined  by 
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Because  of  the  transverse  isotropy  this  also  defines  the 
tranverse  shear  stress  strain  relation  in  any  set  of  transverse 
orthogonal  axes. 

Next  we  consider  the  case  of  an  applied  temperature  rate 
•  • 

<t>  throughout  the  composite  with  all  =  0  in  (2.1)  and  (2.6). 

The  symmetry  is  as  before  and  the  repeating  rectangle  is  again 
subject  to  the  boundary  conditions  (2.12)  but  because  of  the 
transverse  isotropy  the  average  strain  rates  in  the  x^  and  x^  di¬ 
rections  must  be  equal.  Therefore,  from  (2.14) 

n2  ~  n3  =  n  ‘  (2.19) 

The  effective  tangent  thermal  expansion  coefficients  at  the 
current  state  are  defined  by 


22 


=  € 


33 


Tt" 


C11 


At' 


(2.20) 


where  T,  A  denote  transverse  and  axial  respectively.  In  general 
these  expansion  coefficients  are  functions  of  current  stress, 
strain  and  temperature. 

The  basic  boundary  value  problem  for  the  repeating  rectangle 
which  must  be  solved  is  shown  in  figure  5c.  The  input  is  $  =  1. 
The  average  strain  rates  are  determined  from 
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0 


(2.21) 


°22  '33 

utilizing  (2.16),  and  the  constant  axial  strain  rate  is  found 
from 

on  =  0  (2.22) 

utilizing  (2.18).  Because  of  (2.19-20)  these  strain  rates 
are  equal  to  the  tangent  thermal  expansion  coefficients. 

Because  of  the  rate  linearity  of  the  problems  considered, 
loading  (2.11)  can  be  combined  with  heating  to  analyze  by  the 
same  procedures  effective  stress-strain  relations  in  the  pres¬ 
ence  of  change  of  temperature  history. 

ANTIPLANE  STRAIN  -  AXIAL  SHEARING 


We  now  consider  the  axial  shear  loading  (2.6b).  Guided  by 
the  elastic  formulation  of  the  problem  (ref.  14)  we  write  the 
displacement  rates  in  the  form 


U1 

u2  (x) 
u3  (*) 


L 1 2X2  +  c  1 3 X  3  +  Hx2,x3) 


£  °  V 

L12xl 


t13Xl  • 


(2.23) 


It  will  be  shown  further  below  that  c£2  an^  r- J 3  are  the  average 
strain  rates  in  the  fiber  composite  cylinder.  The  only  non¬ 
vanishing  strain  rates  associated  with  (2.23)  are 
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e12  =  e12  +  ***'2 

( 2 . 24 ) 

e13  =  el3  +  ^'3  • 

The  strain  rates  (2.24)  are  equivalent  to  a  principal 
shear  strain  rate 


e 


2 

P 


(2.25) 


whose  direction  n  is  specified  by 
tan  (x2,n)  =  e13/e 12  • 

This  will  result  in  a  principal  shear  stress  rate  in  the  same 
direction 

t  =  2G  (2.26) 

P  t  p 

where  is  the  tangent  modulus.  The  stress  (2.26)  is  equivalent 
to  the  Cartesian  components 

o12  =  cos (x2 ,n) 

(2.27) 

o13  =  Tp  sin (x2 ,n) 

It  follows  that  the  only  nonvanishing  stress  rates  are  the 
ones  given  above. 
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It  may  be  shown  that  this  formulation  results  in  a  well 
posed  plasticity  boundary  value  problem  for  stress  and  strain 
rates  when  the  boundary  conditions  are  of  type  (2.1)  with  (2.6b) 
i.e.  axial  shearing. 

We  consider  the  case  of  axial  shearing  in  the  12  plane.  The 
boundary  conditions  (2.1)  then  have  the  form 


=  o^2n1  (2.28) 

=  0  . 

The  displacement  formulation  (2.23)  then  assumes  the  form 

•  •  • 

u1(x)  =  ££2*2  +  *<x2'x3* 

u2 (x)  =  £J2x1  (2.29) 

u3(x)  =  0  . 

For  the  purpose  of  computation  of  the  function  y,  it  is  neces¬ 
sary  to  establish  its  boundary  conditions  on  the  contour  of  a 
repeating  element.  Figure  6  shows  four  joined  repeating  elements. 
Evidently,  ip  obeys  the  symmetry  condition 

iMx2,x3)  =  iMx2,-x3)  .  (2.30) 


t3(s) 

t2(s) 

T3(S) 
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RT 


It  will  be  shown  that  it  obeys  an  antisymmetry  condition  with 
respect  to  the  x 2  axis.  To  see  this,  consider  the  two  upper 
rectangles.  If  the  direction  of  shear  stress  acting  on  any  one, 
the  right  one  say,  were  reversed,  then  these  two  rectangles  would 
be  in  a  symmetric  state  with  respect  to  the  x axis.  Reversing 
the  direction  of  shear  stress  in  the  right  side  rectangle  to  con¬ 
form  with  the  actual  situation  reverses  the  sign  of  the  field  in 
this  rectangle.  The  conclusion  is  that  the  fields  are  antisym¬ 
metric  with  respect  to  the  x2  axis,  thus 

!|/(x2,x3)  =  -tp(-x2,x3)  .  (2.31) 

It  follows  from  (2.28)  that 

( 0  ,  x3 )  =  0 

and  by  the  same  reasoning 
i)Ma2,x3)  =  0  . 

It  follows  from  (2.27)  that 

^,3(x2,x3)  —  — i)j,3(x2,— x3) 
and  therefore 

'f'  /  3  (*2 ' 0)  —  ^  •  3  (x2,a3^  ^  * 

Thus,  the  rate  boundary  conditions  for  \p  on  the  boundary 
of  a  repeating  element  are  summarized  as 
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ip(0,x3)  =  \p  (a2>x3)  =  0 


t|/,3  (x2,0)  =  ij;,3  (x2,a3)  =  0 


(a) 

(b) 


(2.32) 


The  second  boundary  condition  is  inconvenient  for  finite 
element  applications.  An  alternative  condition  is 


°13*x2'0)  =  °i3^x2'a3^  =  0  (2.33) 

which  is  again  a  consequence  of  the  fact  that  x2  must  be  an  axis 
of  symmetry. 

At  the  fiber  matrix  interfaces,  C^.,,  displacements  and 
tractions  must  be  continuous.  This  implies  from  (2.29)  that 

ipm(ci2)  =  ^f  (C12)  (2.34) 

where  m  and  f  denote  fibers  and  matrix.  It  follows  from  (2.27) 
that  the  only  nonvanishing  traction  rate  component  on  C32  is 

T1  =  °12n2  +  °13n3 

where  n2'n3  are  comPonents  of  the  normal  to  C12.  Thus,  the 
second  continuity  condition  is 


Tlm(C12)  Tlf(C12) 


(2.35) 


To  show  that  e°2  is  the  average  strain  rate,  we  recall  the 
average  strain  rate  theorem  (ref.  14)  which  states  that 
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(2.36) 


f  <V, 


+  u^n^ldS 


where  V  is  the  volume  and  S  the  surface.  Applying  this  to  a 
repeating  element  of  unit  thickness  in  the  x^  direction,  figure  4, 


for  the  strain  we  have 


12 

where  n 


2a 


V  f  <aln 

2a3  •'S  1 


2  +  U2n^)ds 


(2.37) 


n2  ate  the  components  of  the  outward  normal  in  the  1, 


directions . 


It  is  seen  that  n^  vanishes  on  the  enclosing  faces 


while  n2  vanishes  on  the  upper  and  lower  faces  and  on  the  faces 


x3  a3* 


We  introduce  (2.29)  into  (2.37)  recalling  the  bound¬ 
ary  conditions  (2.32a).  Without  loss  of  generality  the  x1  origin 
may  be  taken  at  the  bottom  face  which  makes  U2  vanish  on  this 
face.  As  a  result  of  all  this 

*o  /-a3  -a2-a3 

cn  =  ■— — f  /  a9dx  +/  /  dx„dx-,l  =  r.?.,  .  (2.38) 

12  2a2a3iyn  2  3  JQ  JQ 


dx2dx3 


“12 


In  view  of  (2.21)  the  only  nonvanishing  average  stress  rate 
for  the  lo jding  (2.28)  is 


a12  C12 


(2.39) 


The  relation  between  the  average  strain  and  stress  rates  (2.38-39) 
defines  the  stress  strain  relation  of  the  fiber  composite  in  axial 
shear.  To  obtain  this  relation  numerically,  we  start  from  some 
known  state  of  stress  and  strain  in  the  composite  e.g.  elastic. 

We  impose  a  strain  increment  Ac °2  and  write  all  rate  equations 
as  increments.  This  defines  a  Av  and  stress  and  strain  incre¬ 
ments  Ao12,  Ao13,  A e ^ 2  and  These  are  computed  by  numerical 


solution  of  the  plasticity  problem  for  the  repeating  element 
subject  to  boundary  and  interface  conditions  (2.32-35)  for 
increments.  The  resulting  stress  increment  .  ■  7  is  proportional 
to  Its  average  '^2  defines  the  average  stress  increment 

corresponding  to  an  average  strain  increment  The  procedure 

is  repeated  for  other  strain  increments  and  thus  the  stress-strai 
relation  is  obtained  in  step  wise  fashion. 

The  procedures  described  in  this  section  permit,  in  principl 
analysis  of  stress  and  strain  states  in  the  hexagonal  array  model 
for  any  kind  of  loading  anu  temperature  change  as  specified  by 
(2.1),  (2.3)  by  numerical  solution  of  rate  value  problems  for 

one  repeating  element  for  isothermal  generalized  plane  strain, 
thermal  input  with  no  load, and  axial  shear.  Because  of  rate 
linearity  the  rate  fields  can  be  superposed  at  any  current  state 
and,  therefore,  any  loading  and  temperature  variation  history 
can  be  treated.  Actual  execution  of  a  solution  may  involve  very 
considerable  computer  expenses.  Therefore,  the  procedures  are 
applied  in  this  work  only  to  one  dimensional  average  loadings 
such  as  uniaxial  stress  and  pure  shear.  The  problem  of  deter¬ 
mination  of  effective  stress-strain  relations  for  arbitrary 
loadings  will  be  considered  in  section  3  by  different  methods. 

NUMERICAL  RESULTS 

The  analytical  models  developed  and  described  in  the  pre¬ 
vious  section  were  utilized  in  conjunction  with  an  elasto-plastic 
temperature  dependent  finite  element  analysis  to  predict  one  di¬ 
mensional  stress-strain  responses  for  a  metal  matrix  composite. 
The  loadings  considered  included  tranverse  normal  stresses, 
transverse  shear  stresses,  axial  shear  stresses  and  uniform  tem¬ 
perature  variations.  The  composite  considered  was  a  T-50  graphite 


fiber  in  a  2024-T4  aluminum  alloy.  The  fiber  volume  fraction 
was  30%. 

The  finite  element  code  utilized  was  ANSYS,  which  is  a 
proprietary  analysis  code  developed  by  Swanson  Analysis  Systems 
Inc.  (ref.  16).  A  description  of  the  formulation  of  the  plas¬ 
ticity  relations  utilized  in  ANSYS  can  be  found  in  Appendix  A. 

The  elastic  material  constants  used  in  the  analysis  per¬ 
formed  for  the  fiber  and  matrix  are  shown  in  table  1.  The 
fiber  properties  correspond  to  a  Union  Carbide  T-50  graphite 
fiber.  The  matrix  properties  are  representative  of  a  2024-T4 
aluminum  alloy. 

The  finite  element  model  utilized  for  the  numerical  analysis 
is  shown  in  figure  7.  The  validity  and  accuracy  of  this  model 
was  demonstrated  both  numerically  and  analytically. 

The  method  used  for  assessing  the  accuracy  of  the  finite 
element  model  involved  predicting  elastic  moduli  of  the  uni¬ 
directional  composite.  The  moduli  were  computed  by  applying 
loads  (stresses)  to  the  model  and  finding  the  rssoli.ant  dis¬ 
placements  (strains).  The  boundary  conditions  of  -  -orre*  d 
to  the  symmetry  conditions  described  previously. 

Elastic  moduli  predicted  by  a  finite  element  analysis  were 
compared  with  composite  cylinder  assemblage  results.  The  com¬ 
parisons  are  shown  in  table  2.  As  can  be  seen,  excellent  agree¬ 
ment  between  the  two  models  was  demonstrated.  Additionally,  a 
comparison  of  the  results  predicted  with  the  finite  element  model 
shown  in  figure  7  and  a  considerably  more  refined  finite  element 
model  was  made.  The  refined  model  contained  approximately  twice 
the  number  of  elements  and  nodal  points  shown  in  figure  7.  The 
differences  in  predicted  elastic  moduli  were  of  the  order  of  0.5%, 
thus  again  demonstrating  the  accuracy  of  the  finite  element  mesh. 
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The  prediction  of  stress-strain  relations  for  the  metal  matrix 
composite  under  consideration  requires  that  full  stress-strain  re¬ 
lations  for  the  constituent  materials  be  utilized.  In  the  case  of 
graphite  fibers,  in  the  temperature  ranges  to  bo  considered,  these 
relations  are  linearly  elastic  and  temperature  independent.  Hence, 
only  the  room  temperature  elastic  constants  are  required.  The 
aluminum  matrix  material  is  considerably  different,  however. 

The  aluminum  matrix  is  el asto-plastic  and  temperature  depend¬ 
ent.  The  way  in  which  this  is  modelled  in  AXSYS  depends  upon  the 
material  strain  hardening  rule  desired.  If  isotropic  hardening 
is  desired,  the  stress-strain  response  of  the  material  must  be 
represented  as  linearly  segmented  curves  at  various  temperatures. 
The  representation  used  for  isotropic  hardening  in  the  current 
study  is  shown  in  figure  8- 

ANSYS  requires  a  different  form  of  the  stress-strain  relations 
if  kinematic  hardening  is  desired.  Here  only  two  linear  segments 
are  allowed  for  each  of  the  various  temperatures.  The  kinematic 
hardening  data  used  in  the  current  study  is  shown  in  figure  9. 

In  addition  to  the  stress-strain  response,  thermal  expansion 
coefficients  are  required.  A  curve  representing  the  secant  thermal 
expansion  coefficient  as  a  function  of  temperature  for  the  aluminum 
matrix  is  shown  in  figure  10. 

Utilizing  the  finite  element  model,  boundary  conditions  and 
material  properties  described,  various  one  dimensional  stress-strai; 
responses  were  predicted.  Additionally,  thermal  expansion  coeffi¬ 
cients  were  predicted.  These  results  are  described  here  in  some 
detail  as  they  demonstrate  many  interesting  features  of  metal 
matrix  composites. 


Thermal  Expansion 


The  free  thermal  response  of  the  T-3D  '2024-T4  material  system 
is  depicted  in  figures  11  and  12.  The  responses  utilizing  both 
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isotropic  and  kinematic  hardening  for  the  aluminum  matrix  are 
shown.  The  curves  represent  a  thermal  cycle  from  the  stress- 
free  temperature  to  room  temperature  and  back  to  the  stress-free 
temperature . 

The  states  of  microscopic  stress  and  strain  which  exist  at 
room  temperature  can  involve  considerable  matrix  yielding.  This 
effect  is  clearly  demonstrated  by  comparing  the  decreasing  tem¬ 
perature  response  to  the  increasing  temperature  response.  In 
both  the  isotropic  and  kinematic  hardening  solutions,  the  re¬ 
heat  is  totally  elastic  at  the  lower  end  of  the  temperature 
scale.  The  differences  between  the  cool-down  and  re-heat  curves 
at  these  low  temperatures  corresponds  to  the  matrix  plasticity 
developed  during  cool-down. 

At  higher  temperatures  in  the  re-heat  portion  of  the  curves, 
the  differences  between  the  two  hardening  assumptions  are  apparent. 
The  isotropic  hardening  predictions  remain  elastic  up  to  the 
stress-free  temperature  due  to  the  expansion  of  the  yield  surface 
during  cool-down.  The  kinematic  prediction  involves  additional 
yielding  upon  re-heating  since  this  hardening  assumption  shifts 
the  yield  surface  rather  than  expanding  it.  Thus,  the  re-heat 
caused  the  microscopic  matrix  stresses  to  extend  beyond  the 
shifted  yield  stress  space. 

The  minor  differences  between  the  two  loading  assumptions 
during  the  cooling  cycle  are  due  entirely  to  the  different 
material  representations  required  by  ANSYS  for  the  two  harden¬ 
ing  assumptions. 

The  effects  of  applied  mechanical  loading  on  the  thermal 
expansion  of  this  metal  matrix  system  are  demonstrated  in  fig¬ 
ures  13  and  14.  The  matrix  hardening  assumption  used  for  these 
predictions  was  isotropic.  The  reference  state  of  stress  and 
strain  for  these  predictions  includes  the  effects  of  cooling 
from  the  stress-free  temperature  to  room  temperature  and,  hence, 
matrix  yielding  has  occurred. 


In  figures  13  and  14,  the  effects  of  applied  mechanical 
loading  on  thermal  expansion  are  shown.  The  figures  compare 
free  thermal  expansion  strains  with  thermal  strains  developed 
in  the  presence  of  an  applied,  constant  stress  in  the  x,  direc¬ 
tion.  In  each  of  the  figures,  the  curves  have  been  shifted 
such  that  zero  strains  are  depicted  at  room  temperature. 

Figure  13  demonstrates  the  effects  of  the  mechanical  load¬ 
ing  on  transverse  expansion.  here,  the  mechanical  load  produces 
an  increased  thermal  expansion  in  the  loading  direction  and  a 
decreased  thermal  expansion  perpendicular  to  the  mechanical  load. 
The  same  effect  is  seen  in  figure  14.  The  applied  mechanical 
load  produces  a  decrease  in  the  axial  thermal  expansion. 

One  Dimensional  Mechanical  Loadings 

The  effects  of  processing  stresses  were  graphically  depicted 
in  figures  11  and  12  where  the  thermal  loading  curves  differ 
from  the  thermal  unloading  curves.  This  effect  is  also  clearly 
present  in  the  mechanical  stress-strain  response  of  this  material 
system. 

The  effects  of  processing  on  transverse  extensional  loadings 
are  shown  in  figures  15,  16  and  17.  The  three  curves  in  each  of 
these  figures  represent  transverse  tension  and  compression  with 
processing  stresses  and  transverse  tension  from  a  microscopically 
stress-free  state.  Figure  15  represents  the  stress- st rain  response 
while  figures  16  and  17  depict  Poisson  induced  strains. 

In  figure  15,  the  differences  between  tensile  and  compressive 
loadings  with  processing  stresses  are  attributable  to  the  sign  of 
the  residual  processing  stresses.  The  tensile  loading  produces 
some  micro-level  elastic  unloading  at  low  load  levels  while 
the  compressive  load  produces  only  increased  micro-level  plastic 
loading.  Hence,  the  tensile  stress-strain  response  is  stiffer 
and  has  a  higher  yield  stress.  The  tensile  loading  curve  without 


residual  stress  depicts  the  same  initial  modulus  as  the  tensile 
with  processing  response  but  with  a  much  higher  yield  stress. 

This  is  due  to  large  built  in  stresses  when  processing  effects 
are  included. 

The  Poisson  induced  strains  in  figure  16  show  little  dif¬ 
ference  between  tension  and  compression  with  residual  stresses 
while  the  curve  which  does  not  contain  processing  effects  is 
significantly  different.  This  is  primarily  because  the  tensile 
load  without  processing  is  nearly  elastic. 

The  axial  Poisson  strains  (fig.  17)  show  a  very  curious 
trend.  While  the  solution  without  processing  is  again  nearly 
elastic,  the  tensile  loading  curve  with  processing  has  a  re¬ 
versal  in  slope.  This  is  again  a  function  of  the  sign  of  the 
processing  effects.  The  applied  tensile  load  produces  a  con¬ 
traction  axially.  The  residual  stresses  are  large  and  com¬ 
pressive  in  the  fibers.  Thus,  as  the  applied  load  causes  fur¬ 
ther  matrix  yielding,  the  matrix  modulus  decreases  and  the  fibers 
expand  axially  relieving  some  of  the  residual  compressive  load. 
This  produces  the  effect  shown  in  figure  12.  The  same  rationale 
with  axial  expansion  rather  than  contraction  due  to  the  applied 
load  produces  the  large  strains  under  compressive  loading. 

Thermal  Effects  on  Transverse  Extension 


The  stress-strain  response  of  Gr/Al  can  be  expected  to 
vary  with  temperature  since  the  aluminum  matrix  is  temperature 
dependent.  This  is  clearly  demonstrated  in  figures  18  and  19. 
The  initial  state  for  this  applied  X2  direction  extension  con¬ 
tains  the  residual  processing  stresses  and  therefore  the  effects 
of  processing  and  elevated  temperature  are  present. 

In  figure  18,  the  effect  of  elevated  temperature  is  seen  in 
the  reduced  initial  modulus  at  204°C.  The  effects  of  processing 
combined  with  the  elevated  temperature  are  seen  in  the  increased 
yield  stress  at  204°C,  The  increase  in  yield  stress  is  a  con¬ 
sequence  of  the  reduction  of  processing  induced  stresses  as  the 
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temperature  increases.  Thus,  while  the  response  at  room  tempera¬ 
ture  is  producing  considerable  matrix  yielding  in  addition  to  that 
present  from  processing,  the  response  at  204°C  requires  substantial 
additional  load  to  re-load  the  matrix  to  the  yield  stress  levels. 

The  effects  seen  in  the  Poisson  induced  strains  are  similar. 

At  the  higher  temperature,  the  reversal  in  slope  in  the  axial 
Poisson  strain  is  not  present,  while  the  transverse  Poisson  stra? n 
at  204°C  is  substantially  reduced  in  comparison  to  the  room  tempera¬ 
ture  response. 

Effects  of  Temperature  on  Transverse  Shear  Response 

The  stress-strain  response  for  Gr/Al  for  applied  transverse 
shear  loading  is  depicted  in  figure  20  for  three  values  of  uni¬ 
form  temperature.  These  predictions  were  each  made  from  a  micro¬ 
scopically  stress-free  state  and  as  such  do  not  include  any  pro¬ 
cessing  effects. 

These  data  clearly  show  the  effects  of  decreasing  matrix 
modulus  and  yield  stress  with  increasing  temperature. 

Comparisons  of  the  effects  of  processing  effects  on  ele¬ 
vated  temperature  transverse  shear  response  are  contained  in  fig¬ 
ure  21.  Here,  the  presence  of  residual  micro-level  stresses  are 
seen  to  alter  the  predicted  response  considerably.  The  initial 
yield  stress  has  been  reduced  nearly  to  zero  in  both  the  room 
temperature  and  260°C  responses.  The  initial  elastic  portions  of 
the  curves  are  extremely  limited  when  the  residual  stresses  are 
included  but  do  exhibit  the  same  slope  as  the  solutions  without 
processing  stresses. 

It  is  most  interesting  that  the  processing  stresses  in  con¬ 
junction  with  the  effects  of  elevated  temperature  do  not  promote 
an  increased  yield  strength  at  the  elevated  temperature  as  was 
seen  in  the  transverse  extension  loadings  (fig.  18).  The  reason 
that  this  effect  is  not  seen  in  the  transverse  shear  loading  re¬ 
lates  to  the  elastic  modulus  at  the  elevated  temperature.  In  the 
shear  solution,  the  elevated  temperature  was  260°C  while  the 


extensional  solution  (fig.  18)  elevated  temperature  was  only 
204°C.  Noting  figure  8,  the  elastic  matrix  modulus  is  signifi¬ 
cantly  lower  at  260°C  than  at  204°C.  Thus,  the  material  reaches 
the  yield  strain  (Appendix  A)  at  a  much  lower  stress  value  at 
260  °C . 


Effects  of  Temperature  on  Axial  Shear  Response 

The  axial  shear  stress-strain  relationships  for  this  Gr/ 

A1  system  at  three  different  uniform  temperatures  are  shown  in 
figure  22.  These  predictions  were  made  from  a  microscopically 
stress-free  state  and,  hence,  contain  no  processing  effects. 

The  results  are  typical  of  elevated  temperature  response  of  a 
homogeneous  material  where  the  modulus  and  yield  stress  decrease 
with  elevated  temperature.  The  axial  shear  solutions  were  made 
without  processing  effects  since  this  would  have  required  a 
full  three  dimensional  analysis.  This  complication  is  not  a 
function  of  the  analytical  formulations.  It  relates  to  the 
operation  of  the  analysis  code  used. 

Load  Path  Effects 


In  order  to  determine  the  effects  of  path  dependence  of  the 
material  system  under  consideration,  a  number  of  analyses  were 
made  utilizing  different  paths  to  obtain  identical  final  load 
states.  These  included  combinations  of  load  and  temperature 
as  well  as  combinations  of  bi-axial  loading. 

Comparisons  of  different  load-temperature  paths  are  shown 
in  tables  3  through  5.  Tables  3  and  4  were  generated  using  the 
isotropic  hardening  assumption  and  tensile  and  compressive  trans¬ 
verse  loadings  respectively. 

In  table  3,  the  data  indicate  that  relatively  low  loadings 
produce  less  path  dependence  than  higher  loadings.  This  is, 
of  course,  to  be  expected,  since  the  higher  loadings  produce 
more  plastic  effects.  The  data  presented  also  indicate  that  in 
the  range  of  loadings  considered,  superposition  of  mechanical 
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and  thermal  loadings  would  produce  reasonable  accuracy.  This 
would  be  of  tremendous  significance  in  simplifying  the  analysis 
of  this  material. 

The  data  in  table  4  demonstrate  trends  similar  to  the  table 
3  data.  Here,  as  in  the  higher  tensile  loads  data,  the  axial 
strain  data  show  significant  path  dependence.  The  magnitudes 
of  the  strains  are  very  small,  however,  and  the  differences 
may  be  exaggerated  by  the  numerical  errors  inherent  in  the  solu¬ 
tion  process. 

The  effects  of  path  dependence  with  the  kinematic  hardening 
assumption  were  also  evaluated  for  transverse  tensile  loadings 
and  temperature.  The  comparison  of  the  load  path  effects  is  listed 
in  table  5.  The  differences  here  tend  to  be  larger  than  for  simi¬ 
lar  comparisons  made  for  isotropic  hardening  (table  3) .  The  dif¬ 
ference  of  15%  in  significant  in  that  superposition  of  me¬ 

chanical  and  thermal  effects  is  probably  not  reasonable  for  these 
kinematic  data. 

The  comparisons  made  here  have  shown  significant  dependence 
of  the  composite  material  response  on  the  hardening  assumptions 
for  the  matrix  material  when  mechanical  and  thermal  loadings  are 
considered. 

A  final  evaluation  of  path  dependence  was  made  utilizing 
a  bi-axial  transverse  stress  state  and  isotropic  hardening. 

The  results  of  this  comparison  are  shown  in  table  6. 

The  load  path  effects  seen  in  table  6  are  considerably  smaller 
than  the  combinations  of  temperature  and  load  with  the  exception 
of  the  data.  As  before,  however,  the  magnitude  of  the 
strains  are  very  small  in  comparison  to  the  transverse  strain 
components  and  can,  therefore,  be  expected  to  be  computed  less 
accurately.  Thus,  it  is  apparent  that  mechanical  loading  alone 
does  not  promote  any  significant  path  dependence. 


THREE  DIMENSIONAL  STRESS-STRAIN  RELATIONS. 
ISOTHERMAL  CASE 


In  the  present  part  of  the  work,  the  problem  of  establish¬ 
ment  of  the  three  dimensional  stress-strain  relations  is  considered 
in  a  global  sense  on  the  basis  of  the  transverse  isotropy  of  the 
material.  A  set  of  simplified  effective  stress-strain  relations 
is  established  which  exploits  given  one  dimensional  stress-strain 
relations  obtained  either  numerically  such  as  the  ones  derived  in 
I,  or  experimentally.  Thus,  the  heterogeneous  structure  of  the 
composite  enters  indirectly  via  these  one  dimensional  stress-strain 
relations. 

The  matrix  is  an  elastic  plastic  temperature  dependent  material . 

The  usual  assumption  is  made  that  the  strain  increment  d  . .  is  com- 

1  j 

posed  of  an  elastic  part  de  '  .  .  ,  a  plastic  part  d,“. .  and  a  thermal 
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expansion  increment  de^  .  Thus 
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dF<  0 


dt".  .  = 
i] 


de  .  7  =  x.  ,d4> 

ID  ID 


(3.4) 


where 


Sijkl^’^  -  temperature  dependent  elastic  compliances 
g  =  g  (a,  £,  4>) 

F  =  F(o,  £,<())=  1  -  the  yield  function 

a.  .  (£,  £,  ii)  -  tangent  thermal  expansion  coefficients; 

dependent  on  stress,  strain  and  temperature. 

It  is  seen  that  the  matrix  stress-strain  relations  can  be 
written  in  the  form 


d£  =  R(n,  £,  i+i)d£  +  n(o,  ;)d; 


(3.5) 


The  fibers  are  temperature  dependent  elastic  transversely 
isotropic.  Their  stress-strain  relations  are 
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(3.6) 


d£23  "  d023/2GT 
dt"13  =  d013/2GA 

where 

1  -  fiber  axial  direction 

2,  3  -  transverse  direction 

Efl  -  axial  Young's  modulus 

V-  -  associated  axial  Poisson's  ratio 
A 

Et  -  transverse  Young's  modulus 

v  -  associated  transverse  Poisson's  ratio 

Gt  -  Et/2(1+vt)  -  transverse  shear  modulus 

G„  -  axial  shear  modulus 

-  axial  tangent  thermal  expansion  coefficient 
aT  -  transverse  tangent  thermal  expansion  coefficient 


All  of  the  above  properties  are  temperature  dependent. 

Consider  a  large  cylindrical  specimen  of  a  fiber  composite, 
with  the  fibers  aligned  in  generator  direction  x^,  at  current  con¬ 
stant  temperature  $  and  under  current  load  defined  by  the  boundary 
tractions 


VS)  =  o^n. 


(3.7) 


where  are  space  constant  stresses  and  n ^  the  components  of 
the  outward  unit  normal.  The  average  stresses  in  the  cylinder 
are  then  (ref.  14) 
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13 


=  a?  . 
13 


(3.8) 


However,  the  detailed  states  of  stress  and  strain  as  well  as 
the  average  strains  are  functions  of  the  loading  history  cul¬ 
minating  in  (3.7). 

Suppose  that  there  are  added  a  temperature  increment  d: 
and  boundary  traction  increments 


dT .  (S)  =  do®  .n  . 

1  13  3 


(3.9) 


then 


do .  .  =  da?  .  •  (3.10) 

13  13 

Since  the  stress-strain  relations  (3.5-6)  of  the  phases  and 

their  other  governing  equations  are  linear  in  stress,  strain 

and  temperature  increments,  it  follows  that  the  average  strain 

increments  in  the  cylinder  are  linearly  related  to  the  temperature 

increment  and  to  do?,  thus  to  do. ..  These  relations  involve  co¬ 
il  13 

efficients  which  are  functions  of  the  current  state.  It  will  be 
assumed  that  the  current  state  and  stress  and  strain  enter  through 
their  averages.  Thus 

d£  =  R*(0,  £,  ;)  do  +  r*(o,  T,  ;)  df  .  (3.11) 

It  is  furthermore  assumed  that  the  usual  plasticity  concepts 
apply  to  the  average  states.  Thus,  the  average  strain  increment 
can  be  separated  as  (3.1-4).  It  appears  reasonable  to  require 
that  g  now  become  a  tensor  since  the  material  is  anisotropic. 
Consequently  (3.11)  is  written  in  form  similar  to  (3.5).  Thus 
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Here  S*  .  .  are 

xikl 


while  elastic  unloading  occurs  when  dF  <  0. 
the  temperature  dependent  effective  elastic  compliances  of 
the  composite.  Since  the  composite  is  transversely  isotropic, 
they  are  defined  by  stress-strain  relations  of  type  (3.6)  between 
average  strain  increments  d7^_.  and  average  stress  increments 


do.  .  with  effective  properties,  E*, 
1  J  A 


v*,  etc.  The  yield  function 
F  depends  in  general  on  average  stress,  strain  and  strain  har¬ 
dening  parameters.  The  functions  gijmn  depend  in  general  on 
average  stress  and  strain  history  and  aV j  are  the  temperature 
dependent  effective  tangent  thermal  expansion  coefficients.  In 
the  following,  overbars  denoting  average  stress  and  strain  will 
be  omitted. 

In  order  to  establish  explicit  stress-strain  relations 
a  number  of  simplifying  assumptions  will  be  made: 


(1)  The  average  plastic  strain  t^in  the  fiber  direction 
is  negligible:  The  basis  for  this  assumption  is 
the  large  stiffness  of  the  elastic  fibers  which 
inhibit  plastic  matrix  strain  in  the  fiber  direction 
Thus 
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(2) 


(3) 

(4) 


(5) 

(6) 


The  average  stress  an  the  fiber  direction  does 

not  change  the  state  of  plastic  average  strain: 

The  basis  for  this  assumption  is  the  same  as 
above . 

The  plastic  volume  change  of  t:.e  n.atiix  is 
negligible . 

The  composite  is  initially  t.ansversoly  isotropic 
and  so  remains  during  plastic  flow:  The  second 
part  of  this  assumption  cannot  be  laterally  true 
but  the  alternative  induces  prohibitive  difficulties. 
The  yield  function  F  depends  on  average  stress  and 
on  temperature. 

The  g. .  are  functions  of  average  stress  and  of 
temperature . 


In  the  following  we  shall  establish  a  specific  set  of  stress- 
strain  relations  for  the  isothermal  case.  We  shall  then  proceed  in 
the  next  section  to  discuss  temperature  dependence.  Since  the  fi¬ 
bers  are  elastic,  it  follows  from  (3)  that 
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and  therefore  from  (1) 


^ 22  +  *33  "  ° 


(3.14) 


From  (3.12)  the  average  plastic  strain  increment  is 
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It  follows  from  (3.13-14)  that 
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A  similar 

conclusion  has  been  reached  in  reference 

13. 

The 

specific 

nature 

of  the 

functional  dependence  of  (3. 

21) 

on  the 

invariants  is  not  known.  We  shall  make  the  usual  assumption 
that  F  is  a  quadratic  function  of  average  stress.  This  is 
consistent  with  the  Mises  yield  function  of  isotropic  plasticity 
and  the  Hill  yield  function  of  orthotropic  plasticity.  This 
eliminates  I5  since  it  is  a  cubic  and  leads  to  the  form 


F  =  I3/;T2  +  I4/;a2  =  1 


(3.22) 


where  iT  and  t  are  two  temperature  dependent  constants  with 
dimension  of  stress  which  change  with  the  loading.  It  is  seen 
that  I3  and  1^  appearing  in  (3.22)  are  the  squares  of  the  prin¬ 
cipal  shear  stresses  in  the  transverse  x2x3  plane  and  the  axial 
plane  in  the  direction,  respectively.  Hence,  the  initial 
yield  surface  is 


V’to 


V'ao 


=  1 


(3.23) 


-36- 


where  t_  and  t,  are  the  initial  transverse  and  axial  shear 
To  Ao 

yield  stresses,  respectively. 

If  it  is  further  assumed  that  the  yield  function  grows  in 
stress  space  in  self  similar  fashion,  which  is  the  analogue  of 
the  isotropic  hardening  assumption  for  initially  isotropic 
materials,  then 
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(3.25) 


In  order  to  define  the  stress-strain  relations  it  is  neces¬ 
sary  to  determine  the  forms  of  the  functions  9^jmn  appearing 
in  (3.12).  The  plastic  strain  increments  as  defined  by  these 
equations  are  given  by 


de"  .  .  =  g.  .  -§i- 

lj  ^ijmn  da 


(3.26) 


It  is  recalled  that  the  composite  is  assumed  to  be  trans¬ 
versely  isotropic  throughout  the  process  of  elasto-plastic  defor¬ 
mation.  Since  the  elastic  and  the  free  thermal  parts  of  the 
strain  increment  obey  transverse  isotropic  symmetry  separately 
such  symmetry  must  also  be  valid  for  the  plastic  strain  incre¬ 
ment  (3.26).  Recalling  that  the  yield  function  F  is  not  a  function 
of  Oj^  the  transversely  isotropic  form  of  (3.26)  can  be  written 
as 
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Introducing  (3.28)  into  (3.27)  we  have 
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This  result  can  be  applied  to  the  case  ~  ‘"°33  w^ich  produces 

a  pure  shear  in  a  system  of  axes  rotated  at  45°.  Because  of 
transverse  isotropy,  the  transverse  shear  stress-strain  relation 
for  de^-j  must  aPP^Y  to  this  case. 
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This  procedure  easily  determines  the  stress-strain 
for  ^£23*  Thus  we  have 

dt“22  =  2gT  31^  !5(o22“°33)  dF 

de33  =  ~  2gT  !5(a22"°33)  dF 

de23  =  2gT  7l^  °23  dF 


Denoting 


2g1212  ~  gA 

it  follows  from  transverse  isotropy  that 
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The  plastic  work  increment  dW^  is  given  by 


dW  =  0 .  .d£" .  . 
P  13  1: 


Inserting  (3.31-32)  into  this  expression,  we  find 
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Because  of  the  transverse  isotropy  dWp  can  depend  on 
stress  only  through  the  invariants  of  the  stress  tensor. 

It  follows  that 

gT  =  gT  d3*  i4) 

(3.35) 

gA  =  gA  ^3'  I4* 

To  obtain  information  about  the  functions  gT  and  g  ,  we  shall 
use  one  dimensional  stress-strain  relations  in  transverse  and 
axial  shear.  It  will  be  recalled  that  a  similar  approach  is  use¬ 
ful  in  isotropic  plasticity  for  isotropic  J2  theory.  In  that 
case,  however,  there  is  a  single  function  g  which  is  assumed  to 
depend  only  on  the  single  invariant  .  The  functional  dependence 
g(J2>  is  then  easily  determined  in  terms  of  a  single  stress-strain 
relation  e.g.  uniaxial  stress  or  pure  shear.  Another  variant  of 
this  approach  is  in  terms  of  a  relation  between  so  called  effective 
stress  and  strain  rate  which  are  essentially  square  roots  of  J2  of 
stress  and  strain  rate. 

Hill  (ref.  17)  has  attempted  to  extend  these  concepts  to  or¬ 
thotropic  plasticity  (see  also  further  discussion  in  reference  18) . 
However,  unlike  the  isotropic  case,  the  validity  or  even  the  exis¬ 
tence  of  a  relation  between  such  defined  effective  stress  and  strain 
rate  is  not  evident,  as  has  also  been  pointed  out  by  Hill.  Here 
we  have  to  bring  in  the  anisotropy  by  two  different  functions  gT 
and  gA  or  in  general  by  a  tensorial  g^jmn*  It  should  be  noted, 
however,  that  with  such  a  description  it  is  not  evident  that 
Drucker’s  postulate  will  be  universally  satisfied. 

Returning  to  the  problem  under  consideration,  let  the  prin¬ 
cipal  shear  stress  for  the  plane  stress  system  c22,  °33'  °13  '3e 
denoted  PT  and  the  principal  shear  stress  for  the  stresses 

a.-,  o. .  be  denoted  n  It  is  seen  from  (3.20)  that 
12'  13  A. 
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(a) 

(b) 


(3.36) 


where  superscript  a  denotes  stress  invariant.  We  may  define 
similar  invariants  of  plastic  strain  increment.  Thus 

1^ "  =  ~  d>:33)2  +  <d-23)2  = 

(3.37) 

lJC"  =  (de12)2  +  (dc13)2  =  (bd,A)2 

where  the  extreme  right  sides  are  squares  of  principal  plastic 
shear  strain  rates. 

The  angle  between  pT  and  the  x2  axis  is  given  by 

ta"26T,  '  -  '  '3'38» 

while  the  angle  6^,  between  d  i^,  and  the  x2  axis  is  given  by 

22  -  dc 33 

tan2-TY  =  -  -  He*-  -  (3'39) 


It  is  seen  from  (3.31)  that  these  angles  are  the  same.  There¬ 
fore,  the  directions  of  pT  and  hdy"  coincide.  The  same  conclusion 
is  reached  in  similar  fashion  for  the  directions  of  p  and  ^dr" . 

A  A 

Introducing  (3.31)  into  (3.37)  yields  the  relations 


d'T  *  2«t  fr;  dF 

(3.40) 

dvA  -  2%  Hr  PA  dr 
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With  choice  of  the  yield  function  (3.22)  these  relations 
assume  the  form 


2gT 

<*Y"  =  — j  PT  dF 
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(3.41) 


Consider  the  case  of  isothermal  transverse  principal  shear  pT 
in  the  absence  of  axial  shear.  The  one  dimensional  stress- 
strain  relation  may  be  written  in  the  Ramberg-Osgood  form 


[1  + 


(3.42) 


where  GT  is  the  transverse  shear  modulus  and  sT  and  m  are  curve 
fitting  parameters.  Because  of  the  transverse  isotropy  this  is 
the  same  as  the  stress-strain  relation  for  o23  and  e 23 .  It  fol¬ 
lows  from  (3.42)  that 
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(3.43) 
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Also  in  the  present  case 
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Introducing  (3.4  3b—  44)  into 
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gT(i3,o)  = 


mir 


8gtst 


1-1  P1 


(3.41a) 

-3 


we  have 


Therefore,  in  view  of  (3.36a) 
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(3.45) 


The  last  equation  follows  from  (3.25a)  with  1^  =  0. 

We  now  proceed  in  similar  fashion  for  the  case  of  isothermal 
axial  principal  shear  in  the  absence  of  transverse  shear.  The  Ram- 
berg-Osgood  presentation  for  one  dimensional  axial  shear  is 


,  /  A.  n-1, 

’»  '  S  ”  V  1 


(3.46) 


where  is  axial  shear  modulus  and  and  n  are  curve  fitting 
parameters.  Because  of  the  transverse  isotropy  this  is  the  same 
as  the  stress-strain  relation  for  and  L 12'  We  then  obtain 

by  the  same  procedure  as  before 
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(3.47) 


The  last  equation  follows  from  (3.25b)  with  =  0. 

The  functions  (3.45)  and  (3.47)  define  the  isotropic  har¬ 
dening  plastic  strain  rates  for  plane  stress  and  for  general 
axial  shear  but  not  for  a  combination  of  the  two.  In  the  first 
case,  from  (3.31),  (3.22)  and  (3.45) 
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In  the  second  case,  using  (3.32),  (3.22)  and  (3.47) 

n-1 

dth  -  357  ‘tV  0 12  dF 

A  O. 


(3.49) 
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dc;3  ■  357  ‘75»  =13  dF 


The  preceding  results  are  special  cases.  For  an  arbitrary 
state  of  stress  gT  and  g^  are  functions  of  1^  and  I 4,  (3.35). 

In  attempting  to  construct  such  functions  we  shall  require 
that  in  the  special  cases  discussed  above  these  functions  reduce 
to  (3.45)  and  (3.47).  Furthermore,  stress- strain  relations 
established  must  include  isotropic  J ^  theory  as  a  special  case 
when  the  material  becomes  isotropic  in  shear.  This  will  be  the 
case  when 
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ST  SA  1 To  1 Ao 

Gt  =  Ga  n  =  m  . 

In  this  event 

gT  =  gA 

and  is  a  linear  function  of  and  1^.  To  comply  with  these 
requirements,  a  convenient,  albeit  non-unique,  choice  is  to  assume 
that 


gT  gT  +  ^"4^SA  ^ 

gA  =  gA  (I3/st2  +  I4/sA2)  ' 


(3. 50) 


The  only  functions  (3.50)  which  comply  with  (3.45)  and  (3.47) 
are 
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(3.51) 


Introducing  these  functions  and  the  yield  function  (3.22)  into 
(3.31-32)  and  taking  into  account  (3.22)  we  have  for  the  iso¬ 
thermal  case 
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where 


i ,  j  =  2 ,  3 
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(3.52) 


(3.53) 


S12  ~  °12  S1 3  _  k13 

Equations  (3.52-53)  define  the  incremental  isothermal  plastic 
stress-strain  relations  in  terms  of  the  initial  yield  stresses 
tTo  and  tAo/  the  Ramberg-Osgood  parameters  sT,  m,  sft  and  n  and 
the  elastic  shear  moduli  G  and  G.  . 

i.  n 

We  shall  now  consider  some  special  cases  of  interest- 
1.  Proportional  Loading 

This  implies  that  all  stresses  grow  in  fixed  ratios 
to  one  another.  Thus 
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(3.54) 


where  (o)  is  some  initial  stress  state.  If  (3.54)  are  intro¬ 
duced  into  (3.52-53)  taking  into  account  (3.20),  the  stress-strain 
relations  can  be  integrated.  The  results  are 
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i»  j  =  2,3 

2.  Plane  Stress  in  the  x^x2  Plane 

This  case  is  of  importance  for  balanced  symmetric 
laminates  which  consist  of  unidirectionally  reinforced  layers. 
For  membrane  type  loading,  i.e.  force  resultants  in  the  mid-plane 
of  the  laminate  and  no  bending  moments,  the  state  of  stress 
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in  any  lamina  is  plane  except  for  a  boundary  layer  near  the 
edges.  The  stress  state  on  a  typical  lamina  element  is  shown 
in  figure  23.  It  follows  from  (3.20)  that  in  the  present  case 
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Then  from  (3.52-53) 
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For  the  case  of  proportional  loading  these  equations  integrate 
into 
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Next  we  consider  the  loading  paths  shown  in  figure  24.  The  loading 
OC  represents  the  proportional  loading  considered  above.  The  other 
load  paths  represent  uniaxial  stress  first  and  shear  second,  and 
vice  versa.  For  the  loading  OAC  equations  (3.56)  can  be  integrated 
to  give  the  results 
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For  the  loading  OBC 
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TEMPERATURE  DEPENDENT  STRESS-STRAIN  RELATIONS 


It  is  a  well  known  fact  that  the  yield  stresses  of  materials 
decrease  with  rising  temperature  and,  therefore,  it  is  to  be  ex¬ 
pected  that  yield  surfaces  will  contract  when  the  temperature  in¬ 
creases.  This  implies  that  the  same  state  of  stress  produces 
larger  plastic  strains  for  higher  temperatures.  Phillips,  ref¬ 
erence  19,  has  conducted  extensive  tests  to  determine  the  vari¬ 
ation  of  the  yield  surfaces  of  metals  with  temperature.  His 
results  for  aluminum  indicate  that  the  initial  yield  surface  at 
any  temperature  is  accurately  described  by  a  Mises  ellipse  and 
that  the  dimensions  of  the  ellipse  decrease  linearly  with  rising 
temperature,  figure  25.  We  are  not  aware  of  similar  experiments 
for  metal  matrix  composites  and  in  the  absence  of  such  informa¬ 
tion  we  shall  extend  these  findings  to  the  present  situation  in 
generalized  fashion.  We  shall  assume  that  with  temperature  rise 
the  yield  surface  shrinks  in  self  similar  fashion.  This  implies 
that  all  stress  parameters  defining  the  yield  surface  change  in 
similar  fashion  with  temperature.  In  particular,  the  quadratic 
initial  yield  function  (3.23)  may  be  considered  as  a  yield  func¬ 
tion  at  a  reference  temperature  $  and  is  accordingly  written 

VTTor  +  VTAor  =  1 

Then  the  initial  yield  functions  at  some  other  temperature  is 

13/tTo(*)  +  VTAo(4,)  =  1  (4‘2) 

where  temperature  self  similarity  is  expressed  by 

TTo(4,)  =  TTorf(,p) 

(4.3) 

T  -  ( <t> )  =  T,  f  ( <*  ) 

Ao  T  Aor  Y 
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If  the  yield  stresses  decrease  linearly  with  temperature  as  in  the 
experimental  results  shown  in  figure  26  and  if  the  yield  surface 
"shrinks  to  zero"  at  the  same  temperature  t}>u ,  figure  26,  then 


(4.4) 


At  any  temperature  <p  it  is  assumed  as  in  section  3  that  the 
material  hardens  in  transversely  isotropic  fashion.  Therefore,  the 
subsequent  yield  surfaces  at  temperature  $  are  again  given  by 

I3/tt(<}))2  +  I4/taU)2  =  1  (4.5) 

where  tt  ( <J> )  and  r  (<{>)  are  related  to  tTo  ( )  and  tAq(<J>)  by  (3-25>- 
We  can  also  express  the  temperature  dependence  of  the  initial  yield 
surface  in  the  form 
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Similarly,  the  temperature  dependence  of  subsequent  yield  surfaces 
may  be  given  by 


-  f2  (4>) 
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(4.7) 


where  xTr  and  iAr  are  the  parameters  of  subsequent  yield  surfaces 
at  the  reference  temperature  <(>r .  The  relations  between  tTj_  and 

tt,  ta  is  given  by  (4.3).  Thus 


xT  ( 4> )  =  TTyf(<t>) 
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In  regard  to  the  functions  and  gA  we  can  follow  through  the 
method  of  their  construction,  discussed  in  the  previous  section,  at 
any  temperature  0  with  Ramberg-Osgood  parameters  for  that  temperature. 
Thus,  if  the  shear  stress-strain  relations  at  temperature  0  are  ex¬ 
pressed  as 
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then  in  analogy  to  (3.51) 
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The  stress-strain  relations  are  now 
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and 


where  the  elastic  strain  increments  are  given  by  (3.6), 
are  axial  and  transverse  temperature  dependent  tangent  thermal 
expansion  coefficients,  5^^  is  the  Kronecker  delta  and 
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using  the  yield  function  F  =  0.  Applying  this  to  (4.5)  and  using 
(4.8)  we  find 
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The  temperature  dependent  plastic  strains  in  (4.11)  can  now  be 

written  explicitly  as  follows 
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where  we  have  used  (3.25)  and  (4.3),  and  all  material  parameters  in 
(4.17-18)  are  functions  of  temperature. 

We  consider  some  special  cases  of  (4.13-14). 

1.  Heating  at  constant  stress 
In  this  event 
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When  the  temperature  is  changed  from  $  to  i  the  strains  at 
cf>  are  given  by 
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These  integrals  can  be  carried  out  if  the  temperature  depend¬ 
ence  of  all  physical  parameters  in  (4.9)  is  known.  In  order  to 
determine  such  temperature  dependence,  it  is  necessary  to  have  a 
set  of  transverse  and  axial  shear  stress-strain  relations  at  vari¬ 
ous  temperatures  and  to  model  each  by  the  representation  (4.9). 
Such  stress-strain  relations  can  be  obtained  experimentally  or 
numerically.  Here  we  have  chosen  the  second  alternative.  The 
stress-strain  relations  at  various  temperatures  have  been  obtained 
by  the  methods  described  previously  with  the  help  of  the  ANSYS 
computer  code,  based  on  matrix  property  changes  with  temperature. 
Such  results  were  shown  for  transverse  shear  in  figure  20  and  for 
axial  shear  in  figure  22.  Each  of  the  stress-strain  relations 
was  fitted  into  a  Ramberg-Osgood  form  (4.9)  and  the  resulting 
parameters  sT,  s^,  m  and  n  were  plotted  as  functions  of  tempera¬ 
ture.  Such  plots  are  shown  in  figures  27-30.  The  initial  yield 

stresses  and  t,  were  determined  from  each  stress-strain  curve 

To  Ao 

as  the  points  where  the  linear  relation  ended.  Plots  of  these  as 
functions  of  temperature  are  shown  in  figures  31-32. 

All  of  the  stress-strain  relations  are  based  on  the  premise 
that  the  composite  has  no  internal  stresses  at  room  temperature 
(RT) .  This,  however,  is  not  necessarily  realistic  since  a  metal 
matrix  composite  is  manufactured  at  elevated  temperature  at 
which,  it  may  be  assumed,  there  are  no  significant  internal 
stresses.  As  the  composite  is  cooled  down  to  RT,  significant 
internal  stresses  develop,  tensile  in  the  matrix  and  compressive 
in  the  fibers.  If  the  composite  is  now  loaded  in  shear,  the  re¬ 
sidual  matrix  stresses  interact  with  matrix  stresses  due  to  load 
and,  therefore,  the  stress-strain  relation  will  be  different  from 
that  for  specimens  with  no  residual  stresses  at  RT.  Note  that  the 
linear  part  of  the  stress-strain  relation  is  not  changed  by  the 
residual  stresses. 

To  illustrate  these  phenomena,  we  have  performed  numeri¬ 
cal  analyses  of  the  hexagonal  array  model.  The  composite  has 
been  assumed  stress-free  at  371°C.  The  internal  stresses  due  to 
cool-down  to  RT  have  been  determined.  After  that,  two  subsequent 
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analyses  were  performed.  In  the  first,  a  transverse  shear  load¬ 
ing  was  applied.  In  the  second,  the  composite  was  re-heated  to 
260°C  and  then  the  shear  loading  applied.  Both  these  analyses 
were  shown  in  figure  21  and  compared  with  analyses  where  the 
material  was  presumed  stress-free  at  room  temperature  and  260°C. 
These  curves  demonstrate  that  the  plastic  strains  are  much  larger 
when  residual  stresses  are  included.  In  both  analyses  where 
residual  stresses  were  included,  the  elastic  portion  of  the  re¬ 
sponse  was  severely  limited  in  extent.  The  initial  moduli  were 
identical  to  the  solution  without  processing  stresses,  however. 

Similar  phenomena  can  be  expected  for  axial  shear  but  we 
have  not  performed  the  numerical  analysis  for  this  case. 

These  findings  are  of  significance  for  the  general  stress- 
strain  relations  as  developed  in  this  work.  The  point  of  view 
taken  here  is  that  the  basic  information  entering  into  the  stress- 
strain  relations  is  transverse  and  axial  shear  stress-strain  re¬ 
lations  as  a  function  of  temperature.  If  the  three  dimensional 
stress-strain  relations  are  required  for  a  composite  which  has 
residual  stresses  due  to  cool-down,  then  it  is  necessary  to  use 
the  temperature  dependent  parameters  of  shear  stress-strain  re¬ 
lations  with  residual  stress  in  the  same  fashion  as  done  above. 

This  completes  the  discussion  of  the  three  dimensional  tem¬ 
perature  dependent  stress-strain  relations.  It  is,  of  course, 
necessary  to  examine  the  validity  of  the  stress-strain  relations 
and  this  can  be  done  experimentally  or  numerically.  In  the  lat¬ 
ter  case,  it  is  necessary  to  compute  stress-strain  relations  un¬ 
der  some  states  of  combined  stress  by  the  incremental  methods 
developed  in  the  section  METHODS  OF  ANALYSIS  and  to  compute  the 
same  stress-strain  relations  on  the  basis  of  the  theory  developed 
in  subsequent  cnapters,  based  on  numerically  obtained  one  dimen¬ 
sional  stress-strain  relations,  examining  whether  the  two  sets 
of  stress-strain  relations  are  close.  Unfortunately,  the  scope 
of  the  present  contract  did  not  permit  such  an  examination  since 
the  computer  time  needed  is  quite  extensive.  The  ANSYS  computer 
code,  which  we  have  judged  most  suitable  for  our  purposes,  is 

very  expensive  to  run  for  such  cases. 
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CONCLUSION 


We  have  established  a  general  method  of  analysis  for  incre¬ 
mental  loading  and  heating  of  metal  matrix  composites  modelled 
by  a  periodic  hexagonal  array  of  identical  circular  fibers.  The 
matrix  is  temperature  dependent  isotropic  elasto-plastic  and 
the  fibers  are  temperature  dependent  transversely  isotropic  elas¬ 
tic.  The  analysis  must  be  carried  out  numerically  but  it  has  been 
shown  that  it  is  sufficient  to  analyze  one  two  dimensional  repeat¬ 
ing  element  for  generalized  plane  strain  and  for  antiplane  strain. 
Utilizing  this  method,  we  have  obtained  a  variety  of  one  dimen¬ 
sional  effective  stress-strain  relations  at  various  temperatures 
and  for  various  temperature  histories.  Thermal  expansion  char¬ 
acteristics  without  and  with  external  load  have  also  been  deter¬ 
mined.  In  all  numerical  analyses  we  have  employed  the  ANSYS  code. 

Next  we  have  established  a  general  set  of  incremental  tempera¬ 
ture  dependent  elasto-plastic  stress-strain  relations  for  the  fiber 
composite  based  on  transverse  isotropy  of  the  composite  and  the 
assumption  of  transversely  isotropic  hardening.  A  unique  feature 
of  the  stress-strain  relations  is  that  the  only  information  re¬ 
quired  is  one  dimensional  temperature  dependent  stress-strain  re¬ 
lations  in  transverse  and  axial  shear.  Such  stress-strain  rela¬ 
tions  we  have  obtained  numerically  but  it  is  equally  possible  to 
utilize  experimentally  obtained  stress-strain  relations. 

It  seems  that  the  most  problematic  assumption  in  this  work 
is  transversely  isotropic  hardening  of  the  composite  which  implies 
that  the  yield  surface  changes  in  self  similar  fashion  in  stress- 
temperature  space.  It  is  believed  that  such  an  assumption  is 
sufficiently  accurate  for  initial  loading  of  a  metal  matrix  fiber 
composite  but  not  necessarily  for  repeated  loadings.  Experience 
with  metals  shows  that  for  such  loadings  the  yield  surface  for 
subsequent  loadings  deforms  and  translates  in  stress  space  and 
it  is  possible  that  similar  phenomena  will  occur  for  metal  matrix 
composites. 


Another  idealization  which  is  used  for  metals  is  kinematic 
hardening  which  assumes  that  the  yield  surface  translates  in 
strass  space  without  changing  its  form.  This,  however,  cannot 
be  compatible  with  temperature  dependence  for  the  yield  surface 
must  shrink  with  increasing  temperature  and,  thus,  temperature 
history  cannot  be  accommodated  within  pure  kinematic  hardeninq. 
Thus,  idealized  as  the  assumption  of  transversely  isotropic 
hardening  may  be,  it  is  difficult  to  see  a  viable  alternative 
to  it. 

The  validity  of  the  stress-strain  proposed  must,  of  course, 
be  examined  and  this  should  preferably  be  done  by  experiment. 
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Table  1.  Elastic  Constants  (Room  Temoerature) 


Elastic 

Constant 

T-50  Graphite  Fiber 

2024-T4  Al  Matrix 

E,  (GPa) 

388.2 

72.4 

Et  (GPa) 

!  7.6 

j 

|  72.4 

. 

gat  (GPa) 

14.9 

27.2 

Gtt  (GPa) 

2.6 

27.2 

V 

AT 

0.41 

0.33 

V  TV- 

J  -1 

0.45 

0.33 

r»  m/m/°C 

-0 . 68xl0-6 

22 . 5xl0-6 

uT  m/m/°C 

9.74x10“6 

22.5xl0-6 

A  =  Axial  (Longitudinal) 
T  =  Transverse 
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Table  2.  Elastic  Results  (Room  Temperature) 


Elastic 

Constant 

CCA* 

*  ★ 

F.E. 

% 

Difference 

Et  (GPa) 

41.78 

42.26 

1.1 

Gtt  (GPa) 

14.99 

15.13 

0.9 

G  (GPa) 

A 

22.87 

23.20 

1.6 

VTT 

0.394 

0.396 

0.5 

VAT 

0.338 

0.340  + 

0.6 

aA  m/m/°C 

6.36 

6.36 

0.0 

aT  m/m/°C 

25.65 

25.69 

0.2 

A  =  Axial  (Longitudinal) 

T  =  Transverse 

* 

Composite  Cylinder  Assemblage 

** 

Finite  Element 

+Computed  using  Eft  from  Rule  of  Mixtures 
and  F.E.  Results 
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Table  3.  Transverse  Tension  and  Temperature, 
Isotropic  Hardening 


Load  Path 

£  % 

2 

e  % 

3 

£1  % 

-0.205 

-0.616 

-0.074 

EM 

-0.190 

-0.619 

-0.073 

%  Difference 

7.6 

0.5 

1.4 

0.373 

-0.506 

-0.010 

EH 

0.336 

-0.465 

-0.030 

%  Difference 

10.4 

8.4 

NA 

Load  Steps 

1.  Stress  Free  Temperature  to  Room  Temperature 

2.  Room  Temperature  to  204 °C 

3.  o2  -  0  to  113.8  MPa 

4.  Simultaneous  o2  -  0  to  113.8  and  Room  Temperature  to  204°c 

5.  o2  -  0  to  172.4  MPa 

6.  Room  Temperature  to  299  °C 

7.  Simultaneous  o2  -  0  to  172.4  MPa  and  Room  Temperature 


to  299°C 
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Table  4.  Transverse  Compression  and  Temperature, 
Isotropic  Hardening 


Load  Path 

£2  % 

£3  % 

E1  % 

-0.790 

-0.382 

0.008 

EH 

-0.794 

-0.371 

-0.018 

%  Difference 

0.5 

2.9 

NA 

Load  Steps 

1.  Stress  Free  Temperature  to  Room  Temperature 

2.  Room  Temperature  to  204 °C 

3.  a2  -  0  to  -113.8  MPa 

4.  Simultaneous  02  -  0  to  -113.8  and  Room  Temperature 

to  204°C 
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Table  5.  Transverse  Tension  and  Temperature, 
Kinematic  Hardening 


Load  Path 

e  2  % 

E  ^  % 

E1  % 

-0.201 

-0.614 

! 

-0.079 

1-4 

-0.234 

-0.616 

-0.077 

%  Difference 

i 

15.2 

0.3 

2.6 

Load  Steps 

1.  Stress  Free  Temperature  to  Room  Temperature 

2.  Room  Temperature  to  204 °C 

3.  a2  -  0  to  113.8  MPa 

4.  Simultaneous  -  0  to  113.8  and  Room  Temperature 

to  204 °C 
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Table  6.  o 2  and  a 3,  Isotropic  Hardening 


Load  Path 

e  2  % 

c3  % 

ci  % 

1-2 

0.698 

-0.579 

-0.0214 

2-4 

0.693 

-0.580 

-0.0197 

3 

0.696 

-0.581 

-0.0203 

Max  %  Difference 

0.7% 

0.3% 

7.9% 

Lead  Steps 

1.  Oj  -  0  to  180  MPa 

2.  a3  -  0  to  -90  MPa 

3.  Simultaneous  a2  -  0  to  180  MPa  and  a 3  -  0  to  -90  MPa 
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Figure  2.  Conposite  Cylinder  Assemble: 


omposite  Cy 


ijure  7.  Finite  Element  Model 
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Figure  8.  Input  Matrix  Properties,  Kinematic  Hardening 


Matrix  Properties,  Isotropic  Harde 
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Figure  11.  Transverse  Thermal  Expansion,  Isotropic  vs.  Kinematic  Hardening 
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Figure  13.  Effect  of  Applied  Stress  on  Transverse  Thermal  Expansion,  Isotropic  Harde 
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MPa 


Axial  Strain  vs.  Transverse  Stress, 
Isotropic  Hardening 
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MPa 


Figure  18.  Transverse  Tensile  Stress-Strain  Response  at 
204°C,  Isotropic  Hardening 
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Figure  19.  Transverse  Normal  und  Axial  Strain  vs 
Transverse  Tensile  Stress  at  204°C, 
Isotropic  Hardening 
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Figure  25.  Initial  Temperature  Dependent 

Yield  Surfaces,  Aluminum  (Ref.  19) 
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Figure  27.  Temperature  Dependence  of  Ramberg-Osgood  Parameter 
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Figure  30.  Temperature  Dependence  of  Ramberg-Osgood  Parameter 
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Figure  31.  Temperature  Dependence  of  Initial  Yield  Stress 


of  Initial  Yield  Stress 


APPENDIX  A 
ANSYS  PLASTICITY 

The  finite  element  analysis  code  utilized  in  the  current 
study  was  ANSYS.  This  is  a  proprietary  computer  cod<_  developed 
and  maintained  by  Swanson  Analysis  Systems,  Inc.  ANSYS  is  a 
general  multipurpose  finite  element  code  with  many  capabilities. 

Two  different  hardening  rules  were  used  in  the  metal-matrix 
composite  analysis.  Botli  isotropic  hardening  and  classical  bi¬ 
linear  kinematic  hardening  were  considered  for  the  aluminum  ma¬ 
trix  material.  The  development  of  these  two  hardening  rules 
is  fundamentally  different  in  ANSYS,  even  though  the  yield  con¬ 
ditions  (von  Mises)  and  flow  rules  (Prandtl-Reuss)  are  identi¬ 
cal.  A  brief  synopsis  of  the  pertinent  elements  of  the  deriva¬ 
tions  of  the  two  plasticity  approaches  is  given  here.  A  more 
detailed  description  can  be  found  in  reference  16. 

ISOTROPIC  HARDENING 

In  ANSYS,  isotropic  hardening  in  two  or  more  dimensions  is 
developed  as  an  extension  of  a  one-dimensional  case.  The  pro¬ 
cess  involves  computing  equivalent,  one-dimensional  stresses 
and  strains  for  evaluating  the  yield  condition.  The  procedure 
is  as  follows. 

First,  a  one-dimensional  strain  which  includes  plasticity 
effects  is  computed  at  iteration  i. 


e,i 


T  <l+ve)  | 


l  .  t  t.  ' 
1  x~cy 


,  t  . t. 
( c  - 1  ) 

y  z 


<^x> 


3  2  3  t  2  |  1/2 
+  2('>yz)  +  2(>xz)  \ 


(A-l) 


where  f  =  total  strain-thermal  strain-origin  shift  strain.* 


An  artificial  strain  which  accounts  for  stress  reversal , 
defined  later. 
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This  strain  is  compared  to  the  previous  maximum  emax  and  loading 
is  occurring  if  e  .  >  c  .  If  this  is  the  case,  then  e  is 

©  ,  .1  JnaX  IUua 

set  equal  to  c  •  .  Based  on  e  ,  and  the  input  stress-strain 
e  ,  i  e  ,  l 

curves,  a  one-dimensional  equivalent  stress,  c  .,  is  computed. 

e  f  i 

In  equation  A-l,  vg  is  defined  as: 


v 

e,  i 


1 

2 


e ,  1 


E  c. 


e,i 


(A-2) 


thus , 

e,  l 

iteratively 


and 


c  .  are  related  and  the  computation  must 
6  ;  1 

A  quantity  (d"e/de p) ^  is  now  conlPuted: 


be  made 


doe 

p  l-l 


e ,  i  L  e  ,  i  - 1 

c  .  -  c  , 

e , l  e , l-l 


1 

E 


(*-3) 


which  is  the  rate  of  change  of  equivalent  stress  to  plastic 
strain.  utilizing  this,  the  plastic  strain  increment  can  be 
calculated . 


df 


2  1  +  ■' 

et  ~  3  ~~E  e,i-l 

|d-v)  U+glg^J  .  > 

p ,  i- 1 


(A-4) 


where 


-  4  Ju'-cV  ♦ 


»  »  O  »  •  O  *3  1  0 

-  .'(£-£  )*  +  (c  -c  )*  +  (c  -e  )  +  4<y  ) 

et  3  J  x  y  v  z  zx  2xy 


3  '  ,2  3  ,  '  ,  2  (  1/2 

+  2  1 yz  +  2  ^xz  | 


(A-5) 
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and  e  .  =  e  -  plastic  strain. 

3  3 

With  Ac  p/  eet»  and  cj  defined,  the  plastic  strain  com¬ 
ponents  are  computed  from  the  Prandtl-Reuss  flow  rules. 


it*'  =  ! -  (2c 

x  3cet 


ieP  =  (2c’ 

V  3cet  y 


c  -  c  ) 

y  z 


i  1  \ 

e  -  e  ) 
z  x 


:.cp  =  - 

. -P  _ 

t.  — 

Lc\ 

z 

X 

/.  c 

,  c 

p 

~'xy  " 

3iet 

'  xy 

1  t 

.  p 

P 

-  '  yz 

3t  et 

’yz 

..,P  = 

> 

P 

t 

” '  xz 

3t  „ 

1  XZ 

et 

(A-6) 


If  unloading  is  occurring,  based  on  e  comparisons, 

rnci  x 

then  another  equivalent,  elastic,  strain  is  computed. 


I -  !uf-c*)2  ♦  2  * 


» 2 (1+v) 


\  t-  «- 

I  x  y 


,  e  e, 

( G  — G  ) 
Z  X 


3/e123/e^23/e12-  1/2 
+  2(Yxy>  +  2  Yxz  +  ?  1  y 2  J 


(A-7) 


where  c^  are  elastic  strains. 


The  equivalent  stress  during  unloading  is  simply: 


°e,i  =  E  Es  • 


(A- 8) 


If  emax  <  0,  the  unloading  logic  is  complete.  The  condi¬ 
tions  required  for  the  negative  l  are  described  later. 

max 

Since  the  equivalent  stresses  and  strains  used  in  the 
derivation  are  all  positive,  special  provisions  are  made  for 
evaluating  load  reversal  during  the  unloading  process.  Specif¬ 
ically,  the  largest  elastic  normal  strain  component  is  compared 
to  its  corresponding  plastic  strain  component.  If  they  have 
opposite  signs  then  the  material  is  assumed  to  be  unloading. 

To  account  for  this,  certain  corrections  to  the  strains  are 
made.  These  include: 

«.  =  -  |c  I 

max  I  max  I 

Current  origin  shift  strains  =  previous  origin  shift 
strains  +  2  previous  plastic  strains 

Current  plastic  strains  =  -previous  plastic  strains 


The  effects  of  these  adjustments  is  to  convert  the  unload¬ 
ing  into  loading  by  shifting  the  origin  and  reversing  the  sign 
of  the  plastic  strains.  As  was  stated  previously,  the  isotropic 
hardening  used  is  an  extension  of  one-dimensional  plasticity. 
This  manifests  itself  in  the  stress  reversal  evaluation  and 
therefore  the  isotropic  hardening  is  not  recommended  for  general 
three-dimensional  plasticity.  The  analyses  made  under  the  iso¬ 
tropic  hardening  rules  have  been  shown  to  be  consistent  within 
themselves  and  with  the  kinematic  hardening  analyses,  however. 

KINEMATIC  HARDENING 

The  classical  bi-linear  kinematic  hardening  used  in  ANSYS 
does  not  suffer  the  reversed  loading  difficulties  found  in  the 
isotropic  hardening  behavior  analysis.  This  is  because  it  is 

-104- 


w 


not  developed  as  an  extension  of  a  one-dimensional  case. 

For  the  kinematic  hardening,  the  plastic  strain  increment 
is  defined  differently.  Specifically,  the  e  (eqn.  A -5)  term  is 

•  sh  ■ 

defined  in  terms  of  -  e instead  of  and  the  yield  stress 


is  used  instead  of  o 
sh 


e,i 


sh 

The  strain  c .  is  defined  such  that 
! 


Ee j  equals  the  increase  in  stress  above  the  yield  stress.  In 
addition,  the  term  (doe/dc  )  is  defined  as  E(S)/1-S,  where 

(S ) E  is  the  slope  of  the  second  line  on  part  of  the  input  bi¬ 
linear  stress-strain  curves.  When  the  computed  value  of  Ac 
(eqn.  C-4)  Is  greater  than  zero,  yielding  is  occurring.  Other¬ 
wise,  unloading  is  occurring.  A  specific  example  given  in 
reference  16  for  the  various  strain  calculations  is: 


(ee) .=(£).-  Aep 
x  i  x  i  x 


(cP^)  .  =  (eP^)  .  .  +  Aep 
xi  x  i-l  x 


(A-10) 


Oi-  A 


sh  • 

The  term  Ee..  is  the  amount  that  the  yield  surface  is 
shifted  within  the  stress  space,  in  the  j-direction.  Since  the 
relations  do  not  use  the  equivalent  stresses  or  strains  for 
evaluating  the  plasticity  effects,  reversed  loading  behavior  is 
handled  directly  and  adjustments  required  for  isotropic  harden¬ 
ing  are  not  needed. 
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